All (qubit) decoherences: Complete characterization and physical implementation 
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We investigate decoherence channels that are modelled as a sequence of collisions of a quantum 
system (e.g., a qubit) with particles (e.g., qubits) of the environment. We show that collisions induce 
decoherence when a bi-partite interaction between the system qubit and an environment (reservoir) 
qubit is described by the controlled- (7 unitary transformation (gate). We characterize decoherence 
channels and in the case of a qubit we specify the most general decoherence channel and derive 
a corresponding master equation. Finally, we analyze entanglement that is generated during the 
process of decoherence between the system and its environment. 



I. INTRODUCTION 

One of the most distinctive features of quantum sys- 
tems is their ability to "exist" in superpositions of mutu- 
ally exclusive (orthogonal) states [1]. Providing a quan- 
tum system has been prepared in a pure state then 
we can write = J2k c k\ t l ) h)i where \ipk) are orthonor- 
mal vectors that compose a basis ((V'ilV'fe) = $ki)- All 
bases are unitarily equivalent and we can express the 
same state in different bases. In fact, we can always 
select a basis such that |^) is a basis vector so in its 
matrix representation the vector is represented by a 
single diagonal element. According to quantum postu- 
lates for the isolated system any evolution is governed 
by unitary transformations and the original information 
about the state preparation of the quantum system is pre- 
served. As soon as an interaction with an environment 
comes into the play (the quantum system is open) the 
situation becomes dramatically different and the state 
is no longer described by the single diagonal element in 
some basis. Depending on properties of the environment 
and the character of the interaction our system evolves 
non-unitarily and its state is, in general, described by 
a statistical mixture. Among various possible dynamics 
of an open quantum system interacting with its environ- 
ment a specific role is played by a process in which the 
off-diagonal elements of the original state g = |vE , )(^E'| in 
some basis are continuously suppressed in time, i.e. 



Q Qt- 



diag[e]. 



(1.1) 



This is a process of decoherence during which some of the 
information about the initial state of the quantum sys- 
tem might be irreversibly lost [2-4]. The basis in which 
the decoherence takes place is specified by properties of 



the environment and the character of the interaction [4] . 
There are at least two aspects of quantum decoherence 
that keep it in the center of interests in multiple investi- 
gation related to foundations of quantum mechanics and 
in quantum information processing. The first aspect is, 
that decoherence is presently viewed as a mechanism via 
which classicality emerges from the realm of quantum 
(see e.g. [2-6]). In this context it is of paramount im- 
portance to specify the basis (the so called pointer basis 
[4]) in which the decoherence takes place. In the field of 
quantum information the decoherence is an evil - it de- 
grades quantum resources (superpositions of states and 
quantum entanglement) that are needed for quantum in- 
formation processing [7] . The degradation of resources is 
caused by random interactions (errors) between a quan- 
tum system under consideration (e.g. a qubit or a quan- 
tum register) with its environment. If nothing else then 
these two facets of quantum decoherence are enough to 
justify an investigation of decoherence channels (trans- 
formations). 

As mentioned above the decoherence is caused by (un- 
avoidable) interactions between the system and its en- 
vironment. Consequently, the whole process of decoher- 
ence can be completely described within the framework 
of the quantum theory as a unitary process that gov- 
erns the joint evolution of the quantum systems and its 
environment 1 [2-4] . There are plentiful theoretical mod- 
els describing the decoherence within the framework of 
the standard quantum theory that have been in accor- 
dance with various experiments [9,10]. These models ei- 
ther use Hamiltonian evolution of the composite system- 
plus-enviroment structure (the Hamiltonian itself is time- 
independent). 

Alternatively, the description of decoherence can be 
based on a simple collision-like models, i.e. a sequence of 



'Another possibility would be to include decoherence into the basic dynamical equation, i.e. to add a non-hamiltonian part 
into the Schrodinger equation [8]. However, the modifications of the basic quantum dynamical law are out of scope of this 
paper. 
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interactions between the object under consideration and 
particles from environment leads to decohcrence. These 
models allow us to study microscopic dynamics of open 
systems, in which the flow of information from the sys- 
tem to the environment and creation of entanglement can 
be analyzed. In fact, collision models are equivalent to 
more general models of causal memory channels [11]. In 
this case, the memory is represented by the system under 
decoherence, whereas the reservoir (environment) plays 
the role of input/output systems. 

In the present paper we will focus our attention on 
collision-like models of decoherence of qubits. Our first 
aim is to completely classify all possible decoherence 
channels of a qubit. The second task is to show that 
all decoherence maps of qubits can be modelled as se- 
quences of collisions. The paper is organized as follows: 
Sections II and III are devoted to a description of gen- 
eral properties of all decoherence channels. In Sec. IV 
we present a generic collision-like model. In the Sec. V 
the master equations for collision models are derived and 
all possible master equations describing decoherence of a 
qubit are presented. In Sec. VI we analyze how entangle- 
ment is created during a sequence of collisions. Finally, 
in Sec. VII we summarize our results and formulate some 
open problems. 

II. DECOHERENCE CHANNELS 

The aim of this section is to classify all possible com- 
pletely positive trace-preserving maps (quantum chan- 
nels) that describe quantum decoherence. Let us denote 
by V the set all maps £ satisfying the decoherence con- 
ditions, i.e. 

(e k \£[g]\e k ) = (e k \g\e k ) for all k (2.1) 
|(e fe |£[f>]|e ; >| < |<e fc |e|e,)| for all k^l, (2.2) 

with B = {\e k )} being the decoherence basis. For our 
purposes it is useful to fix one basis B and to analyze all 
decoherences (forming the set V B ) with respect to this 
basis. The general decoherence maps are then just uni- 
tary rotations of elements from V B , that correspond to 
a change of the decoherence basis. In particular, if £ is 
a decoherence map, then also £' = IA\£U2 is such a map. 
We used the notation Uj[g] = UjgUj with Uj unitary op- 
erators. ^From the definition it is clear that decoherence 
channels are unital (they preserve the total mixture, i.e. 
£[I] = I) and are not strictly contractive (they might 
have more than a single fixed point). 

Denoting by T> B the set of all decoherence maps with 
respect to a fixed basis B we can write T> = U B T> B . Each 
decoherence map £ € T> belongs only to one class T> B . 
Elements of T> B and T> B > arc unitarily related, i.e. 

V B , = {£'\£'[g] :=£[UgU\£ EV B ,B' = UB} = V UB . 

This defines a new decoherence class only if B' ^ B. 
That is, the unitary operation U does not commute with 



all projectors \e k ){e k \, or equivalently the basis B is not 
an eigenbasis of the transformation U. If [U, |efc)(efc|] = 
for all k then from a given £ £ T> B we obtain different 
decoherence maps within the fixed set V B . 



A. Qubit decoherences 

In what follows we will analyze the case of qubit deco- 
herence channels. In this case the set T> has surprisingly 
simple form. We will use the so-called left-right nota- 
tion, in which the evolution map is represented by a 4x4 
matrix [12]. Let us choose the following operator basis 

5 = I; 

51 = IVXV^I + I^XV'I ; 
s 2 = i|V)(V ± |-# ± )(V'l ; 

53 = 1^)^1-1^X^1, (2-3) 

where B = {\"4>}, \ip ± )} is the decoherence basis. The 
elements of 5-basis satisfy the same properties as the 
Pauli operators, because Sj — Wo-jW^ with W be- 
ing a unitary operation. In this basis the operators 
(states) take the form of four-dimensional vectors g = 
\{I + r-S)<r+r e = {l,r), where rj = Tr[gSj]. The evolu- 
tion £ is described by 4x4 matrix with elements given by 
the equation £ k \ = ^Tr(S k £[Si\). Because of the trace- 
preservation we have £oo = 1 an d £qi = £02 = £03 = 0. 
Consequently, we obtain the Bloch sphere representation 
[7] of the state space, in which the states are illustrated 
as points (three-dimensional real vectors r) lying inside 
a sphere with a unit radius. The action of £ corresponds 
to an affine transformation of the Bloch vector r, i.e. 
r — > f 1 = Tr + t, where Tj k = £j k (for j, k = 1,2, 3) and 
tj = £jo. The translation vector t describing the shift 
of the Bloch sphere (including its center, i.e. the total 
mixture) is related to the unitality of the channel. For 
unital maps t = 0. 

Diagonal elements of the state g are in this case associ- 
ated with the mean value z = Tr[gSs\. The conservation 
of the diagonal elements implies that the corresponding 
components of g are preserved. Combining the unitality 
with this property we find the following form for deco- 
herence maps 

/ 1 \ 

«- l: b A • («> 

\o 1/ 

from where it follows that the set of all possible qubit 
decoherence maps is at most four-parametric. 
Each unital map can be written as [12] 

£[Q] = Ru^sRuM - U^ £ [U 2 gUl]u\ , (2.5) 
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where Ru 1 , Rjj 2 are orthogonal rotations corre- 
sponding to unitary transformations U\,U 2 ] = 
diag{l, Ai, A2, A3} and A^ are the singular values of the 
matrix £ . In fact, the above relation is the singular-value 
decomposition of the matrix £. The conditions of the 
complete positivity restricts the possible values of Xj. 
In particular, the allowed points A = (Ai,A2,A3) must 
lie inside a tetrahedron with vertices that have coordi- 
nates (1,1,1), (1,-1,-1), (-1,1,-1), and (-1,-1,1), 
respectively. 

Applying these facts to the decoherence map under 
consideration (£ from Eq.(2.4)) we obtain that $g = 
diagjl, Ai, A2, 1}, i.e. A3 = 1. Let us note that in this 
case we use unitaries that do not change the decoherence 
basis, so we are still dealing with all decoherences that 
belong to a fixed basis B. The condition of complete pos- 
itivity restricts the values to the points A = (A, A, 1) with 
— 1 < A < 1, i.e. to a line connecting the two vertices of 
the tetrahedron representing the identity (A = 1) and the 
unitary rotation S3 (A = —1). Consequently, the general 
decoherence channel £ € T>b reads 



/ 1 

d si 

-si ci 

Vo 1 



/ 1 

A 

A 

\ 1 




(2.6) 



where Sj = sin ipj and Cj — cos tpj represent rotations 
Rjj. around the z-axis by an angle <fj. ^From here it 
follows that a general decoherence map £ takes the form 



£ = 



10 

Xcos(ipi + ip 2 ) \sm(ipi + (p 2 ) 

-Asin(<^i + ip 2 ) Acos(<^i + ip 2 ) 

1 



(2.7) 



and, consequently, it is specified only by two real param- 
eters a = A cos(iy9i + ip 2) and b = A sin(<^i + ip 2 ), i.e. 



£ = 



( 1 
a b 
-b a 

Vo 1 



(2.8) 



As a result we obtain that any map £ of the above form 
with the numbers a, b satisfying the condition a 2 + b 2 < 1 
is completely positive. Therefore we can conclude that 
the set of all decoherence maps of a qubit is characterized 
just by two parameters. Moreover, to obtain the deco- 
herence (to secure the suppression of off-diagonal terms) 
the inequality must be strict, i.e. a 2 + b 2 < 1. Otherwise 
the map £ desribes a unitary rotation around the z axis. 
Defining the rotation map 



cos f sin if 
— sin if cos f 



(2.9) 



and using the relation <p = fi+f> 2 , we can write the most 
general decoherence channel (£ € T>b) in a very compact 
form 



/ 1 

\ R 

^ 

V 1 



(2.10) 



This form is suitable for our purposes, because the pow- 
ers of the map £ read 



£ n 




(2.11) 



III. STRUCTURAL PROPERTIES OF 
DECOHERENCE CHANNELS 

In this section we will briefly review structural prop- 
erties of the set of all possible decoherence completely 
positive maps £. Let us denote this set by V. 

• Convex structure 

The set of all decoherence maps T> is not convex, 
i.e. a convex combination of two decoherence chan- 
nels £ fl = [i£\ + (1 — n)£ 2 is not again a decoherence 
channel. This is true except the case when the de- 
coherence bases of £i,£ 2 coincide, i.e. the set T>b is 
convex. The extremal points of T>b correspond to 
unitary transformations. However, these are not el- 
ements of T>b , because they do not fulfill the second 
decoherence condition (2.2). 




FIG. 1. (Color online) The cube corresponds to all positive 
unital trace-preserving maps. The condition of complete pos- 
itivity confines quantum channels into the tetrahedron with 
(generalized) Pauli matrices as vertices. In this picture the 
set of decoherence channels T>b forms a line connecting the 
points / and S z . 
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We have already mentioned that for qubits the set 
of all possible $g channels form a tetrahedron and 
up to unitary trasnformations each channel belongs 
to this tetrahedron. Those channels that corre- 
spond to dccoherence maps form a line connecting 
the points (1, 1, 1) and (—1, — 1, 1). From this pic- 
ture (see Fig. 1) the convexity of T>b is transparent 
and also the extremal points can be easily identified 
as unitary channels. It follows that each decohcr- 
ence map can be written as a convex sum of only 
two unitary channels. In fact all maps for which 
one of the A's equals to unity and all the others are 
the same define a decoherence with respect to some 
basis. This means that all edges of the tetrahedron 
correspond to decoherence channels. It illustrates 
that the set V as a whole is not convex, but is com- 
posed of a continuous number of "convex" subsets 
T>s corresponding to each orthonormal basis B. 

• Composition 

A composition of two decoherence channels £ = 
£\ o£ 2 is not, in general, a decoherence channel. So 
the set V is not closed under the operation of mul- 
tiplication. The channel £ belongs to V only if the 
decoherence bases of £\ and £ 2 coincide, i.e. again 
only the sets T>b are closed under the composition. 

• Classical capacity 

The decoherence basis is preserved by the decoher- 
ence map. Therefore it is possible to exploit these 
bases states to transmit the maximally possible 
amount of information, i.e. the capacity achieves 
its maximum C = log 2 d with d = dimH. 

• Tensor product 

The tensor product of two decoherence maps de- 
scribes a decoherence. However, T>i2 ^ V\ <g) 2? 2 , 
because the decoherence basis of £\ <S> £ 2 is always 
separable. The open problem is whether the whole 
set £>i2 can be obtained from the sets Vi,V 2 by 
global unitary rotations. Properties of decoherence 
channels under tensor products is an interesting 
topic, which is related to our ability of controlling 
the decoherence. For example, how the decoher- 
ence of a sub-system affects characteristics of the 
whole system? 



IV. COLLISION MODEL 

In what follows we will study whether an arbitrary de- 
coherence channel can be implemented via a sequence of 
bi-partite collisions. Each of the collisions is described 
by a unitary transformation U. Our task will be to de- 
rive all possible unitary transformations that force the 
system to decohere. Our analysis will be performed only 
for qubits, but up to technical details all results hold for 
qudits. 



Let us consider that initially the system qubit is de- 
coupled from an environment (reservoir) that is modelled 
as a set of qubits, i.e. fli n = g ® S res . Moreover, we will 
simplify the model by assuming that initially the reser- 
voir qubits are in a factorized state E res — £® N and each 
reservoir qubit interacts with the system qubit just once. 
In addition we assume that reservoir qubits do not in- 
teract between themselves. Under such conditions the 
evolution of the system qubit is induced by the sequence 
of maps £\ = . . . = £jq = £. In particular, the state of 
the system after the n-th interaction equals to 



=£„... £i[g] = £ n [g] 



(4.1) 



where £[g] = Tr res {U(g®£,)W]. We will refer to this pic- 
ture as to to a collision model. The system qubit collides 
with reservoir qubits. 

In order to obtain the decoherence channel, i.e. 



(n) 

£00 012 
(n) 
021 £11 



with g^ = [f>2i^]* ^ for n goes to infinity, we have 
to ensure that the map £ preserves diagonal elements of 
each state g in a given (decoherence) basis. 

In order to preserve the diagonal elements of pure 
states |0)(0| and |1)(1| (decoherence basis) the bi-partite 
unitary transformation U must necessarily satisfy the re- 
lations 



100} -> |(ty> ; 
101} - |0^) ; 
|10)-|1^); 
|ll)-|l</>>. 



(4.2) 



In what follows we will prove our main result that the 
class of possible bi-partite interactions that induce deco- 
herence in collision models coincides with the set of all 
controlled- U transformations (the so-called U-processors 
as introduced in Ref. [13]), where the system under con- 
sideration plays the role of the control and the reservoir 
particle is a target. Certainly, we have to identify those 
transformations for which the off-diagonal elements of 
the system density operator do vanish in the limit of in- 
finitely many collisions with reservoir particles. 

The unitary bi-partite transformation (the controllcd- 
U operation) defined by the relations (4.2) can be rewrit- 
ten into the following operator form 



U= |0)<0|®V + |l)(l|®Vi 



(4.3) 



where Vo,Vi are unitary rotations of a reservoir qubit. 
In particular, V = \ip)(0\ + |V'- L }(1| and V 1 = ^(O] + 
\<fi)(l\- Thus, the initial state = g <g> £ of a bi-partite 
system evolves according to a transformation 



n' = untf = ]T g jk \j) (k\ ® VjtvZ (4.4) 



j,k=0 
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and by performing the partial trace over the reservoir 
qubit we obtain the induced map 

1 

Q^q' = S[q] = Tr p O' = ]T g jk Tr[V^}\j)(k\ 

j,k=0 

= dmg[g] + Qoi(X) € \0)( 1| + eio <X + )*|l)<0| , 

where X = V?V and (A) e = Tr[A£] stands for the mean 
value of the operator X in the state £. 

Applying the transformation £ in a sequence of n col- 
lisions the state of the system qubit is described by the 
density operator 

= = diag[g] + em (X)^|0)(1| + g w (X^\l)(0\ , 

from where we can conclude, that providing |(A}^| < 1 
and < 1 the off-diagonal terms vanish. However, 

because XX* = X^X = /, i.e. X is unitary, its eigenval- 
ues are just complex square roots of the unity. Therefore, 
for the eigenvectors of X the off-diagonal terms do not 
tend to zero. 

The fact that for convex combinations of the eigenvec- 
tors the off-diagonal elements still vanish might sound 
counterintuitive. But it can be seen from the following 
consideration: Let us denote by e IV and e 111 the eigen- 
values of X associated with the eigenvectors and 
I/2), respectively. Then the mean value (X)^ for the con- 
vex combination £ = a|/i)(/i| + (1 — a)|/2)(/2| equals to 
(X)t = e^a+e"'(l-a). The condition \{X}(\ < 1 can be 
rewritten as the inequality 2a(l — a)[l — cos(ip — 77)] < 0, 
which is satisfied only if cos(p — t]) ^ 1, or a ^ 1 and 
a ^ 0. The latter property means that £ is the eigenstate. 
The first property requires ip = t], i.e. the operator X is 
proportional to the identity X — e %ip I . However, under 
this assumption V\ = e llf Vo, i.e. we have no interaction 
and?7 = (e^|0)(0| + |l)(l|)®Vi. Hence, we can conclude 
that whenever the reservoir state is not an eigenstate of X 
and the interaction is not trivial, the described collision 
model with controlled- U interaction forces the system to 
decohere. 

It is straightforward to show that unitary interactions 
U = |0}(0|®y + |l)(l|®V r i induce maps of the left-right 
form (see Eq.(2.8)) with the parameters 



a= 1 -((X) i + (X*) i ); 
b= l -({X)s-(X^), 



(4.5) 



or, equivalently, (X)^ — \e %Lp . So given a decoherencc 
map £ one can, in principle, find an interaction U and 
an initial state of the reservoir qubits £, such that the de- 
sired decoherence process is implemented via a sequence 
of collisions. 



V. MASTER EQUATION 

In this section we will derive a master equation that 
describe the decoherence process induced by collisions of 



the system qubit with reservoir particles. Although the 
studied decoherence model is intrinsically discrete, we 
will show that we can perform a continuous-time approx- 
imation that enable us to write down the master equation 
(see, e.g. [14]). 

As shown in the previous section the collision model 
is described by a set of maps £ n = £ n that form a dis- 
crete semigroup, i.e. £ n £ m = £ n +m for all integer m,n 
and £0 = X. The question is whether we can introduce a 
continuous one-parametric set of transformations £ t such 
that £ tn = £ n for t n = nr (r is a time scale roughly corre- 
sponding to the time interval between two interactions). 
It turns out that a simple relation n — ► t/r can be used 
to accomplish the task. The obtained continuous set of 
transformations £t will be used to derive the generator Q 
of the dynamics by using a simple formula Qt = £t£f ■ 

With the help of results from Sec. Ill (namely, 
Eq.(2.11)) we can directly write 



£t = 



( 1 




V 









\ 




1 J 



(5.1) 



where for simplicity we set the time scale r = 1 . It is easy 
to see that the one-parametric set of transformations £ t 
possesses the semigroup property, i.e. £t£ s = £t+s- for 
all real t, s. It means that the generator and the associ- 
ated master equation will be of the Lindblad form [15], 
i.e. the process under consideration is Markovian. 
The corresponding generator reads 



/ 

In A -ip 

ip In A 

\ 



(5.2) 



where we used the identity R tlfi — pRt v +-K/2 and 
■^[X 1 Rttp]^ 1 R-t v = hiAi?o + This step can be 

performed only if A is non- negative (i.e., when the loga- 
rithm is defined), which, in general, is not the case. The 
parameter A belongs to the open interval (—1, 1). Conse- 
quently, it seems that the generator cannot be derived in 
all cases. However, using the equality —\R V — Ai?( v+7r ) 
for A nonnegative, one can write |A| t i? t ( l)9+7r ) instead of 
\Rtip in the expression for £ t with A < 0. Then the gen- 
erator is slightly different and contains the term ip + tt 
instead of ip, and ln|A| instead of In A. Thought this is 
not a problem, because in terms of parameters of the col- 
lision model (X)% — Xe tv , i.e. the parameter A = 
is always positive. Therefore we can consider the gener- 
ator Q as the most general one. 

The general master equation in Lindblad form reads 



Qt 



G[Qt] = ~i[H, g t ] + - ^2 c ab ([S a , Q t S b ] + [S a Q t , Sb] ) . 



a, b 
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If the numbers c a b are time-independent and form a pos- 
itive matrix, then the generated evolution is Markovian 
and satisfies the semigroup property. To find the values 
of the coefficients c a b we will use the following relations 
(see Ref. [14]) 

, _ [G\z2 ~ [gjgg . , _ [g]i 3 -[g] 3 i . h _ [g]21 - [G]l2 . 
m — 1 ! n 2 — 4 , "3 — 1 j 



e23 



[Gh 



e3i 



[g]ao . 
4 ' 



ei2 



[fife 



(5.3) 



Let us now address the following question: Is there 
any other master equation describing a decoherence of 
a qubit? The preservation of the S z component (deter- 
mined by the decoherence basis) together with the uni- 
tality of the transformation implies that 



/0 

a b 

c d 

\0 



(5.7) 



and 



dn = 



d 2 2 = 



^33 = 



[G]n 


- [G]22 


- [GU 


[GU 


4 

- [G]n 


- [G} 33 


[Ghz 


4 

- [G]xi 




4 



di2 = 



[g]l2 + [g]21 
4 



*23 



_ [g]23 + [g] 3 



[g]l3 + [g] 3 l 

4 



(5.4) 

where [G]u correspond to matrix elements of the gener- 
ator Q, c ab = d ab - ie ab and H = J2 a h a S a - Note that 
d a b form a symmetric matrix and e a b is an antisymmetric 
matrix. 

Using these expressions one finds that the non- 
vanishing parameters are 



l 



In A 



and the corresponding master equation reads 
Qt = -i^[S3, Qt] ~ ^-(#3 ft S3 - Qt) ■ 



(5.5) 



(5.6) 



The corresponding matrix C 



[Cab_ 



then reads 



C 



a — d c + b 
c+b d-a 
-a-d 



(5.8) 



This matrix is positive only when a — d and b = — c. 
Moreover, a must be negative. These restrictions leave 
only a single element that does not vanish, namely, 
C33 = —a/2. Consequently, the Hamiltonian part takes 
non- vanishing value for /13 = 6/2. Therefore the family 
of all master equations describing the decoherence is only 
two-parametric 



Qt = -^[^3, Qt] - 7;( s 3QtS3 - Qt) 



(5.9) 



This general master equation is of the same form as the 
one derived for the collision model (5.6). The parame- 
ters A, tp are related to the parameters of the underlying 
unitary interaction via the formula (X)^ = Xe llp . Let us 
note the constraint A = E [0, 1], since X is unitary. 

Therefore In A < as it is required by the condition on 
possible values of a. 



A typical evolution driven by this equation is depicted in 
Fig. 2. 



VI. ENTANGLEMENT IN DECOHERENCE VIA 
COLLISIONS 



Bloch sphere 




FIG. 2. (Color online) The decoherence of a qubit gov- 
erned by Eq.(5.6). The Bloch sphere that represents the ini- 
tial state space of a qubit is mapped into the line connecting 
the decoherence basis states. On the right the evolution of 
Bloch-vector components for two different initial states is de- 
picted. 



We start with definitions of entanglement quantities 
that we will evaluate. Let us denote the joint state of the 
system of TV + 1 qubits (the system qubit and N reser- 
voir qubits) by £1. The bipartite entanglement shared 
between a pair of qubits j and k can be quantified in 
terms of the concurrence [16] 



C 



max{0, Ai - A 2 - A 3 - A4} , 



(6.1) 



where Xj are decreasingly ordered square roots of the 
eigenvalues of the matrix Rjf~ = QjkOy ® a yQ*jk <J y ® a y 
and Qjk — Tr-j:Q is the state of two qubits under consid- 
eration. 

The case of multi-partite entanglement is a more com- 
plex phenomenon and there is no unique way of its 
quantification. Fortunately, for pure multi-qubit systems 
there is an accepted method of characterization (identifi- 
cantion) of intrinsic multi-partite entanglement. Specifi- 
cally, let us consider how strongly the j-th qubit is corre- 
lated with the rest of qubits in the multi-partite system. 
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This degree of entanglement can be quantified via the 
so-called tangle (see Ref. [17]) 

t 3 = 4dct Qj = 2(1 - Trg 2 ) , (6.2) 

where Qj = Tr-fi is the state of the j-th qubit. Then we 
evaluate bi-partite concurrences between the given j-th 
qubit and any other qubit in the system, i.e. we evaluate 
N quantities Cjk- 

Wootters and his coworkers have found (see Ref. [17]) 
that for pure three-qubit states the inequalities 

5][Cy 2 <T fc Vfc = l,2,3, (6.3) 

hold. In addition they have conjectured that such in- 
equalities also hold for any number of qubits. This con- 
jecture (to so-called Coffman-Kundu- Wootters (CKW) 
inequality) has been recently proved by Osborne [18] 
These inequalities quantify the property which is known 
as the monogamy of entanglement (the entanglement can- 
not be shared freely in multipartite systems). 

As a consequence of the CKW inequality one can de- 
fine a measure of intrinsic multipartite entanglement Aj 
as 

A J = - . (6-4) 

k±j 

where we have used the notation Tjk = [Cjk] ■ It 
is important to note that in the multi-partite case (in 
particular for more than three qubits) the differences 
Aj := Tfc — J2j^k T i k ta ^ c different values for different 
j. Therefore, a weighted sum A = ^X^/A? ^ s an a P~ 
propriate measure of an intrinsic multipartite entangle- 
ment. Based on this quantity we can argue that there 
are multi-partite entangled states for which the entan- 
glement has purely bi-partite origin, as for example the 
family of W states [20] that saturate the CKW inequali- 
ties, i.e. A = 0. 

Let us assume that the system qubit is initially pre- 
pared in the state \x) = a\0) + b\l) and each qubit of the 
reservoir is in a pure state \ip), i.e. the joint initial state 
is |fio) = |x) ® \ip)® N . After n collisions governed by 
bi-partite controlled unitary operations (4.3) the whole 
system evolves into the state 

|fi n ) = [ o |0> ® \V ^)® n + b\l) <g> \V^ n ] ® \^ N ~^ . (6. 

In order to be able to evaluate the entanglement quan- 
tities we have to specify all two-qubits and single-qubit 
density operatos. In particular, for k < n, j < k the 
bi-partite states are given by expressions 

g ok (n) = |a| 2 |0Vo}(0^o| + |&| 2 |l^i)(lV>i| 

+a6*|(7/. |Vi)| ( "" 1) |OVo)(l^i| +cc. ; (6.6) 
Qjk (n) = |a| 2 |V>oV>o)<^o| + \b\ 2 \^i){^i\ ' (6.7) 

where we used the notation |^o) = Vo\ip) and \ipi) = 
Vi\tp)- The single qubit states are as follows: 



eo(n) = |a| 2 |0}(0| + |6| 2 |1)(1| + a6*|(^|^i)r|0)<l| + c.c. 
describes the system qubit after n-th collision, and 

^(n) = |a| 2 |^o)(Vo| + |&| 2 |^i>(V'i| (6.8) 

describes the fc-th qubit of the reservoir after the colli- 
sion with the system qubit. Evaluation of the tangles is 
straightforward and results in expressions 

r (n) =4|a| 2 |6| 2 (l-|(^ |V'i)| 2n ); (6.9) 

r k (n) =4|a| 2 |6| 2 |^o|^ 1 L }| 2 ; (6.10) 

r ok (n) = 4| a | 2 |6| 2 |^o|^i}| 2( "- 1) |(V'o|V'i L )| 2 ; (6.H) 

r J - fc (n)=0. (6.12) 

One can directly verify the validity of the CKW inequal- 
ities 

JV 

Y, r jk (n) = r j0 (n) = 4M 2 KV^i)| 2( "- 1) I(V<o|Vt}| 2 

k=0,k^j 

<4| a 6| 2 |^ |^i L }| 2 -r J (n); (6.13) 

N 

X> fc (n) = nx 4M 2 |(^o|V'i)| 2(n - 1) l(^r)| 2 
fc=i 

<4|a6| 2 (l-|^ |^)| 2n ) = r (n) , (6.14) 

where we have used the relations |(V'o|V'i)l 5= 1 an( I 



i 




2 4 6 8 10 12 14 n 



FIG. 3. (Color online) The behavior of entanglement as a 
function of number n of collisions between the system qubit 
and reservoir qubits. The degree of entanglement between 
•the system qubit and n reservoir qubits after n collisions is 
given by to - it increases with the number of collisions (time) 
to a steady-state value. On the contrary, all reservoir qubits 
after their interaction with the system qubit are entangled 
with the constant degree of entanglement (see the tangle r k ). 
The bi-partite entanglement To k (the square of the concur- 
rence Cofc) is zero until the fc-th reservoir qubit collides with 
the system qubit. After the collision the entanglement takes 
a non-zero value, though it decreases due to subsequent col- 
lisions of the system qubit with other reservoir qubits. It is 
interesting to note that all Tofc(n) for n > k are described by 
the same function. We assune the following initial state of 
the system qubit \ip) = -75 (|0) + |1)) and | (ip Q \tpi) | 2 = 0.75. 
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In the limit of large number of interactions [n — > oo) 
all two-qubit correlations vanish (i.e. finally there is 
no bi-partite entanglement between qubits in reservoir), 
but the entanglement between the system qubit and the 
whole reservoir converges to a finite value 

T ^4|a6| 2 ; (6.15) 
Tbfc -» . (6.16) 

It means that after the process of decoherence the sys- 
tem qubit is not entangled with the reservoir via bipar- 
tite entanglements, but is entangled to the reservoir via 
multi-partite correlations. The final state belongs to the 
family of Grccnbcger-Horn-Zeilinger states that exhibit 
purely multi-partite correlations. 

^From the above one can see how the entangle- 
ment is related to the decoherence. Given the relation 
K^olV'i)! = K-^Ovl ~ ^ wc conclude that the decoherence 
rate restricts the maximum amount of created entangle- 
ment and simultaneously it determines the decrease of 
entanglement with the number of collisions. 

VII. SUMMARY AND CONCLUSIONS 

In this paper we have studied qubit decoherence chan- 
nels as defined by Eq. (1.1). We have presented their 
complete classification. In addition, we have shown that 
all decoherence channels can be modelled as collisions of 
a quantum system with its environment. The bi-partite 
collisions between the system and reservoir particles are 
modelled as the controlled-f/ operations such that the 
system particle is a control while a reservoir particle is 
a target. Using the collision model we have derived the 
most general decoherence master equation in the Lind- 
blad form that describes decoherence. The specific basis 
in which the decoherence takes place as well as the deco- 
herence rates are specified by properties of the controllcd- 
U operation and the initial state of reservoir particles. 
We have shown that in the collision model the decoher- 
ence is accompanied (or, from a different point of view, 
one can say that the decoherence is due to) quantum en- 
tanglement that is created between the system particle 
and the reservoir particles. We have derived the explicit 
expressions for entanglement measures (the concurrence 
between an arbitrary pair of particles involved in the dy- 
namics and a tangle that characterize a degree of entan- 
glement between the given particle and the rest of the 
system). Using these measures and the Coffman-Kundu- 
Wootters inequalities we have shown that in the case of 
decohering qubit collisions between this qubit and the 
reservoir lead to intrinsic multi-qubit entanglement of all 
qubits involved in the process. 

We conclude our paper with some remarks. 

i) 

Even though through the paper we have been paying at- 
tention mostly to decoherence of qubits many of the re- 
sults hold in general. In particular, within the framework 



of a collision model with the controlled-?/ bi-partite col- 
lisions (the system particle plays the role of the control 
while particles from the reservoir are targets) a decoher- 
ence of qudits can be described as well. 

ii) 

The collision model used in this paper is a discrete one. 
We have assumed that a collision between two particles 
is localized in time, so that at a given time instant the 
controllcd-C/ operation (a bi-partite gate) is applied. The 
sequence of interactions is then labelled by an integer 
number n and the total dynamics is represented by a 
discrete semigroup. 

As shown in the paper it is straightforward to intro- 
duce a continuous time parameter so that the continu- 
ous evolution version of the sequence of collision is de- 
scribed by a Markovian process represented by a con- 
tinuous semigroup. We have derived the corresponding 
master equation that describes the process of decoher- 
ence. More importantly, we have shown that for qubits 
this master equation describing the decoherence is unique 
and takes the form (5.9) that can be written as 

g t = -i[H,Q t ]-^[H,[H,g t )}, (7.1) 
2j 

where we use the notation H = |6"3 = and 7 = 
—b 2 /2a — (p 2 /2\n\. We note that the double commuta- 
tor term is well known and usually appears in decoher- 
ence models even for higher-dimensional systems. For 
example, Milburn in his work on intrinsic decoherence 
(see Ref. [8]) has derived a generalization of the usual 
Schrodinger equation exactly in the form (7.1). 
in) 

We have shown that the decoherence in the collision 
model is accompanied (caused) by a creation of entan- 
glement between the system and the reservoir. Unlike 
in the process of homogenization described in [19-21], in 
which the created entanglement saturates the CKW in- 
equalities, in the case of decoherence the entanglement 
results in the Greenberger-Horn-Zcilinger type of corre- 
lations [22]. This means that decoherence process (as 
described by our collision model) does not create an en- 
tanglement between the environment particles. Specifi- 
cally, if we trace over the system qubit (which decoheres) 
in the nth step of the evolution (see Eq.(6.5)), we find 
that the environment is in a separable state 

u en v(n) = Tr|f2 rl }(f2„ 

= [|a| 2 (|Vo)<Vo|) 8 " + |6| 2 (|Vi}<V'i|) 8 "] 

® |V)(VI® (JV " n) , 

where all the parameters are specified in Sec. VI. The 
decoherence rate A and the rotation parameter ip can be 
adjusted by a suitable choice of the interaction U and 
the state of the reservoir £. The collision model reflects 
microscopic origins of both these parameters that enter 
the decoherence master equations. The eigenvalues of 
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the Hamiltonian H are given by the value of ip and the 
parameter 7 is specified by both these parameters. The 
eigenvectors of H form the decoherence basis. 

iv) 

We have shown explicitly that an arbitrary decoherence 
channel for a qubit can be represented via the collision 
model with a particularly chosen controlled-C/ interac- 
tion. However, this result holds for arbitrary dimen- 
sion (i.e. for qudits) as well. Let remind us that an 
arbitary quantum map £ can be represented as unitary 
operation on some larger system (this is a content of the 
Stinespring-Kraus dilation theorem [7]). We have shown 
that for decoherence channels the collision (represented 
by a unitary transformation) must be of the form of the 
controlled-?/ operation. An open question is whether 
each decoherence master equation (even for dim7i = 00) 
can be derived from the collision model. Knowing a de- 
coherence master equation (i.e., knowing a generator Q) 
it is easy to "fix" a time step t = t and define £ T = £. 
This map is for sure a decoherence channel and can be 
realized by a collision U. By applying this "elementary" 
map many times (a sequence of collisions) we obtain a 
discrete semigroup of the powers of £. The inverse task 
is trickier, that is, how do we interpolate between these 
discrete sequence of transformations (parameterized by 
number of collisions) to obtain a continuously parame- 
terized channel. From a construction of the problem we 
know that the solution exists (we have started our anal- 
ysis from the master equation). The question is whether 
this interpolation for qudit channels can be performed as 
easily as for qubits, i.e. by replacing the discrete powers 
of n with continuous parameter t. Nevertheless, given 
the fact that we have started with a continuous set of 
channels £t and by replacing i-trwe obtained £\ = £. 
Consequently, it is possible to replace mi/r to obtain 
the original continuous semigroup of decoherence chan- 
nels £t- As a result we have found that a collision model 
can be used not only to describe any decoherence master 
equation, but can also be used to describe any quantum 
evolution governed by the Lindblad equation. On the 
other hand, it has to be stressed that collision models 
describe evolutions that might not be "interpolated" by 
continuous semigroup of quantum channels 2 . 
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